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Abstract. Convection is an important phenomenon in the atmospheres of A-type and cooler 
stars. A description of convection in ATLAS models is presented, together with details of 
how it is specified in model calculations. The effects of changing the treatment of convection 
on model structures and how this affects observable quantities are discussed. The role of 
microturbulence is examined, and its link to velocity fields within the atmosphere. Far from 
being free parameters, mixing-length and microturbulence should be constrained in model 
calculations. 
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1. Introduction 

The gross properties of a star, such as broad- 
band colours and flux distributions, are signif- 
icantly affected by the effects of convection in 
stars. Consequently, our modelling of convec- 
tion in stellar atmosphere models can signifi- 
cantly alter our interpretation of observed phe- 
nomena. 

1.1. Mixing-Length Theory 

Convection in stellar atmospheres is usu- 
ally based on mixing- length theory (MLT) of 
Boh m-Vitense I lll958h . In this model a sin- 
gle bubble of gas rises a certain mixing-length 
(l/H) before dispersing (Fig. [^). The prob- 
lems with this theory is that it is clearly too 
simple and that the mixing-length is a to- 
tally free paramete r. In their discussion of 
the ATLAS6 models lRelvea & Kuruczl il97«l> 
found discrepancies between theoretical and 
observed uvby colours which might be the re- 
sult of inappropriate treatment of convection 
within the models. Subsequently, several at- 




(a) (b) (c) 



Fig. 1. Schematic bubble representations of convec- 
tion treatments. In mixing-length theory (a), a single 
bubble rises within the atmosphere, while in turbu- 
lent convection bubbles of varying sizes rise (b). In 
(c) we have overshooting above the convection zone. 



tempts h ave been undertake n to improve the 
situation. Lester et al. I dl98 2). for example, in- 
troduced "horizontally averaged opacity" and 
a "variable mixing length" which improved the 
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match with observed uvby colours, but did not 
remove all the discrepancies. 

1.2. The Canuto & Mazzitelli (CM) 
Model 

ICanuto & Mazzitelli I dl991l 1 19921) proposed a 
turbulent model of convection in order to over- 
come one of the most basic short-comings of 
MLT, namely that a single convective element 
(or "bubble" or "eddy") responsible for the 
transport of all the energy due to convection. 
This new model accounts for eddies of vari- 
ous sizes that interact with each other (Fig.^i). 
The CM convection model w as impl emented in 
the ATLAS9 code bv lKupkal i 19961) and has no 
user adjustable free parameters. An improved 
variant is the s elf-co nsistent (CGM) method of 
ICanutoet all i ll 9961) . 

1.3. Convective Overshooting 

Convective bubbles rise above the convections 
zone into the stable regions (Fig. [IJ). This 
is called overshooting, and should be present 
in our model atmosphere calculations. The 
ATLAS9 models introduced an "approximate 
overshooting" whic h has not b een without its 
critics (see Cas te 111 et al. Il997l for full details). 
The following quote from Kurucz' web site is 
aptly summarizes the situation: 

convective models use an overshooting 
approximation that moves flux higher 
in the atmosphere above the top of the 
nominal convection zone. Many people 
do not like this approximation and want 
a pure unphysical mixing-length con- 
vection instead of an impure unphysi- 
cal mixing-length convection 

1.4. Atmospheric Structure 

At T ef{ = 8000K, CM gives essentially ra- 
diative temperature gradient with significantly 
less convective flux than MLT, while ap- 
proximate overshoo ting introduce s flux into 
higher layers (Heit er et al. 1 120021) . Figure [2] 
shows the situation for slightly cooler models 
(r e ff=7000 K). The CM model remains close to 
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Fig. 2. Models with T cff = 7000 K and logg = 4.0 
showing differences arising from changing the treat- 
ment of convection. The upper panel shows the vari- 
ation of temperature with optical depth, while the 
lower panel shows the convective flux as function of 
optical depth 

the radiative temperature gradient. MLT gives 
more convective flux than CM, even when l/H 
= 0.5. Overshooting produces an excess of 
convective flux in higher layers, which pro- 
duces a noticeable bump in the temperature- 
depth relation compared to MLT without over- 
shooting. 

2. The ATLAS CONVECTION Card 

A single control card specifies how the treat- 
ment of convection within ATLAS. 

If the card is not present in the ATLAS 
input, convection is turned off. A couple of 
variables have default values which should be 
noted: mixing-length MIXLTH = 1 and amount 
of overshooting 0VERWT = 1. 
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CONVECTION OFF 

Ensures that convection is turned off for 
the model calculation. Note that this sets 
MIXLTH to 1, but leaves OVERWT un- 
changed. 

CONVECTION ON <MIXLTH> 

Turns on convection within the model 
calculation. The value of mixing-length 
MIXLTH is set to the user specified value, 
but OVERWT is not modified. Hence OVERWT 
will be the default value of 1 and approxi- 
mate overshooting is enabled by default. 
CONVECTION OVER <MIXLTH> 
<OVERWT> 

This variant of the control card, which al- 
lows full control over app roximate con- 
vectiv e overshooting. See Castelli et al. I 
( 1997) for details of convective overshoot- 
ing. 

Note that setting OVERWT = turns off ap- 
proximate overshooting. 
CONVECTION OVER <MIXLTH> 
<OVERWT> <NCONV> 
ATLAS puts the constraint that the convec- 
tive flux (/conv) must be zero above layer 
NRHOX/2 (i.e. higher than the middle layer 
in the model atmosphere desk) This was 
originally introduced to remove numerical 
artifacts. The value is usually a good num- 
ber, except for coolest models, where con- 
straint ge nerates a jump i n T(tr oss ) for T e g 
< 4000K jCastellill20(m In order to alle- 
viate this problem, NCONV introduced into 
a version of ATLAS by Castelli to allow 
user to specify the layer above which / conv 
is surely zero. The default value is 36. This 
is NRHOX/2 for NRHOX = 72 which is typical 
number of layers used in a model calcula- 
tion. 



2.1 . The CM and CGM routines 

Drop-in replacements for the CONVEC and 
TCORR subro utines were implemented by 
lKu^ Jl996l). 

Usage is similar to standard ATLAS 
version: CONVECTION ON <MIXLTH>. 
However, there are differences depending on 
whether the CM or CGM routines are used: 



CM In this case MIXLTH not actually used, but 
must be set to > in order for the routine to 
work correctly. This is to ensure that rest of 
the ATLAS codes knows that convection is 
running. 

CGM Here MIXLTH has the meaning of a* 
(see ICanutoetal.lll996l for details). The 
standard value used is 0.09. 

3. Testing convection models 

None of the current Id models of convec- 
tion are totally satisfactory. 2d and 3d nu- 
merical simula tions are producin g impressive 
results (Frevtag & Steffenl I2004)). as are im - 
proved analytical Id treatments dKupka 120041) . 
However, in order to be confident that the cur- 
rent generation of ATLAS models are produc- 
ing reliable results, we need to know how good 
the treatment of convection is in ATLAS and 
the limitations. 

3.1. uvby photometry 

Small ey" fe Kupkal lll997l) compared the pre- 
dicted uvby colours for the CM model with 
that from the standard ATLAS9 MLT models 
with and without "approximate overshooting". 
Comparison against fundamental T e ff and log g 
stars revealed that the CM models gave bet- 
ter agreement than MLT without overshooting. 
Models with overshooting were clearly dis- 
crepant. This result was further supported by 
stars with obtained from the Infrared Flux 
Method (IRFM) and \ogg from stellar evo- 
lutionary models. However, some discrepan- 
cies still remained, including a "bump" around 
6500 K in the logg obtained for the Hyades 
and continued problems with the Stromgen 
mo index . Similar conclusions were found by 
Schmidt (1999) using Geneva photometry. 

3.2. The mo index 

The mo index is sensitive to metallicity and 
microturbulence, but also convection efficiency 
llRelyea&Kuruczl Il978t l&nalley & Kupka 
119971) . Inefficient convection (CM and MLT 
l/H ~ 0.5) clearly works in the domain of the 
A stars down to b - y ~ 0.20 (-7000 K, F0). 
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Fig. 3. The variation of m index with b-y show- 
ing the sensitivity to microturbulence (£ tul b), mixing- 
length (I/H) and "approximate overshooting". At 
each temperature the model with logg = 4.0, £ tU rb 
= 2 kms~' and l/H = 0.5 is denoted by a square. 
The arrows indicate the effect of varying £ Ulr t,, l/H 
and including o vershooting (for l/H = 1.25). The 
Philip & Esret 1 1980) main-sequence is included 
for reference. 



For cooler stars, convection becomes more ef- 
ficient and substantive within the atmosphere 
and higher values of mixing-length and lower 
microturbulent velocities would be required to 
fit the observed m<j indices of the main se- 
quence (Fig. QJ- The 2d nume rical rad iation 
hydrod ynamics calculations of Ludwi g et al. I 
( 1999) indicate a rise in mixing-length from 
l/H -1.3 at 7000 K to l/H = 1.6 for the 
Sun (5777 K), while a much lower l/H ~ 0.5 
was found at 8000 K ( lFrevtaglfT995T) . This 
is in agreement with that implied by mo in- 
dex. However, around 6000 K there still re- 
mains a significant discrepancy, which could 
only be reduced by invoking the approximate 
overshooting option (Fig. [3}. None of the con- 
vection models used in classical model atmo- 
spheres allows for the reproduction of the mo 
index, unless l/H and the amount of "approx- 
imate overshooting" are varied over the H-R 
Diagram. 

3.3. Stellar Fluxes 

The stellar flux distribution is influenced by the 
effects of convection on the atmospheric struc- 
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Fig. 4. Fluxes for CM models and MLT models, 
with l/H = 0.5 and 1.25, compared to that for a 
model with zero convection. Their are noticeable 
differences, especially in the region 4000 ~ 5000A, 
and the effect of overshooting is considerable. 



ture of the star. As we have seen with photo- 
metric colours, these effects have a clearly ob- 
servable signature. Hence, high-precision stel- 
lar flux measurements should provide signifi- 
ca nt and useful i n form ation on convection. 

iLester et al. I dl982l) presented a study of 
convective model stellar atmospheres using 
a modified mixing-length theory. They found 
small, systematic differences in the optical 
fluxes. Their figures also demonstrate that con- 
vection can have a measurable effect on stellar 
fluxes. 

Figure@]shows the effect of changing mix- 
ing length from 0, through 0.5 to 1.25 on 
the emergent flux for solar-composition mod- 
els with r eff =7000 K, log,g= 4.0 and £ui- b = 2 
kms 1 . The differences are noticeable, with the 
effect of overshooting being considerable. 

3.4. Balmer line profiles 

The temperature sensitivity of Balmer 
lines makes them an excellent diagnos- 
tic tool for late A-type stars and cooler 
llGardiner I 120001). However, as e mphasised 
bv Ivan't Veer & Megessierl dl996l) . the Ha 
and H/3 profiles behave differently due to 
convection: Ha is significantly less sensitive 
to mixing-length than H/3. Both profiles are, 
nevertheless, affected by the presence of 
overshooting, with H/3 being more influenced 
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Fig. 5. The effects of convection on the predicted 
shape of Balmer profiles for models with T c ff = 
7000, logg = 4.0, [M/H] = 0.0 and £ urb = 2 kms~'. 
Ha (upper panel) is unaffected by the values of l/H, 
but sensitive to "approximate overshooting", while 
Hfi (lower panel) is sensitive to both. 



than Ha (Fig. |5}- Since Ha is formed higher 
in the atmosphere than H/3, Balmer lines 
profiles are a very good depth probe of stellar 
atmospheres. Naturally, Balmer profiles are 
also affected by microturbulence, metallicity 
and, for the hott er stars, surface gravity 
llHeiter et al. 120021) . 

In their studi es of Ha and Hjg pro files 
of A and F stars iGardiner et al. I (Il999l) and 
Smallev et al. (2002) found good agreement 
with fundamental stars for CM and MLT 
(l/H ~ 0.5) without approximate overshooting. 
However. iGardiner et al. I i 19991) found l/H = 
1.25 gave better results for cooler stars (T e ff > 
7000 K). 
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Fig. 6. The variation of microturbulen ce with ef- 
fective temperature. Based on results of iGrav et al. I 
1200 ll) for stars near the main sequence. The dashed 
line indicates the approximate variation with T^. 
Note the apparent relatively abrupt change in be- 
haviour between 6500 and 7000 K 



4. Microturbulence 

Microturbulence is a free parameter introduced 
to allow abundances from weak and strong 
lines to agree. It is an extra source of broad- 
ening, which is added to thermal broadening 
of stellar lines. Physically, it is postulated as 
small-scale turbulent motions within the atmo- 
sphere, where the size of the turbulent elements 
is less than the unit optical depth ( lGravll l992). 

Microturbulence does appear to vary 
with effective temperature. Several studies 
have found that microturbulence appe ars to 
vary systematical ly with 7V,ff (Chaffee 1970; 
Nissenl 1198 ll ICouprv & Burkhartl [1992 
Gravetal. 11200 ll) . Fig. [6] shows the variation 



of £turb with T e ff for near main-sequence stars 
(logg > 4.0) bas ed on the results given by 
IGrav et al~l ll200ll) . Microturbulence increases 
with increasing T e g, peaking around mid-A 
type, before falling away to zero for B-type 
stars. There is a relatively abrupt change in 
behaviour between 6500 and 7000 K, which 
is related to the change from weak subsurface 
convection to the fully convective atmospheres 
of cooler stars. 

4.1. Line Asymmetries 

Velocity fields are present in stellar at- 
mospheres which can be measured using 
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Fig. 7. The Convection Recipe for stars near the main sequence ISmallev 1 2004) 



line bisectors (iDravinsI Il987t iGravl 119921) . 
Compared to solar-type stars, the line bisectors 
in A-type stars are reversed, indicating small 
rising colu mns of hot gas an d larger cooler 
downdrafts (Landstreet 1998). It is these mo- 
tions that are thought to be responsible, at 
least in part, for the existence of microturbu- 
lence. In fact, 3d numerical simulations of so- 
lar granulation can account for observed line 
profiles witho ut the need fo r any microtur- 
bulence ( Asplu nd et al. l l2000). Si milar results 
have been found for Proc yon (J^rav 1985; 
lAllende Prieto et al 2002), which is also a 
star with well-known physical parameters (e.g. 
iKervella et al. 120041) . 

Numerical simulations avoid the need 
for a micro tu rbulen ce free parameter 
( Aspl und et al. I 120001) . The microturbu- 
lence of Id is not turbulent motions, but rather 
velocity gradients within the atmosphere. 
Hence, microturbulence should no longer a 
free parameter, but ought to be constrained 
within ATLAS model calculations. Infact, 
Kurucz presented an empirical method for 
constraining depth-dependent microturbuluce 
within ATLAS at this Workshop llKurucz I 
120051) . 



5. A Convection Recipe 

Smallev (2004) presented a schematic varia- 
tion of microturbulence and mixing length with 
r e ff for stars near the main sequence (Fig. [7J. 
For stars hotter than AO there is no convection 
or significant microturbulence. For the early 
A-type stars there is essentially no convec- 
tion within the atmosphere, since the tempera- 
ture gradient is radiative, but there are velocity 



fields as indicated by the modest microturbu- 
lence values. Velocity fields increase as we go 
through mid to late A-type stars, and inefficient 
convection is required within the atmosphere. 
Once convection becomes efficient (F-type and 
later) the value of microturbulence is found to 
drop, while the mixing-length increases. 

6. Conclusions 

The effects of convection on the stellar atmo- 
spheric structure can be successfully probed 
using a variety of observational diagnostics. 
The combination of photometric colours and 
Balmer-line profiles has given us a valuable 
insight into the nature of convection in stars. 
High quality observations that are currently 
available and those that will be in the near fu- 
ture, will enable further refinements in our the- 
oretical models of convection and turbulence 
in stellar atmospheres. 

Overshooting is still an issue to be re- 
solved, since there are clearly velocity fields 
above the convection zone. While the "approx- 
imate overshooting" of Kurucz appears to have 
been discounted by observational tests, there is 
clearly the need for some sort of overshooting 
to be incorporated within model atmosphere 
calculations. 
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